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ABSTRACT 

The structure of magnetic fields within protostellar disks may be studied via polarimetry provided 
that grains are ahgned in respect to magnetic field within the disks. We explore alignment of dust 
grains by radiative torque in T Tauri disks and provide predictions for polarized emission for disks 
' viewed at different wavelengths and viewing angles. We show that the alignment is especially efficient 

in the outer parts of the disks. In the presence of magnetic field, these aligned grains produce 
polarized emission in infrared wavelengths. We consider a simple model of an accretion disk and 
provide predictions for polarization that are available to the present-day instruments that do not 
, resolve the disks and will be available to future instruments that will resolve the disks. As the surface 

■ magnetic field and the bulk magnetic field play different roles for the disk dynamics, we consider 
separately the contributions that arise from the surface areas of the disk and its interior. We find that 

, the polarized emission drops for wavelengths shorter than ~ IQ^Lvn. Between ^ lO/im and ^ 100/im, 

^ ' the polarized emission is dominated by the emission from the surface layer of the disks and the degree 

\ of polarization can be as large as ~ 10% for unresolved disks. We find that the degree of polarization 

OO . at these wavelengths is very sensitive to the size distribution of dust grains in the disk surface layer, 

04 ' which should allow testing the predicted grain size distributions. The degree of polarization in the 

far- infrared/sub- millimeter wavelengths is sensitive to the size distribution of dust grains in the disk 
interior. When we take a Mathis-Rumpl-Nordsieck-type distribution with maximum grain size of 

■ 500-1000 /im, the degree of polarization is around 2-3 % level at wavelengths larger than ~100/im. 
Our study indicates that multifrequency infrared polarimetric studies of protostellar disks can provide 

. good insights into the details of their magnetic structure. 

Oh' Subject headings: accretion, accretion disks — circumstellar matter — polarization — stars: pre-main- 

Q . sequence — dust, extinction 

■ 1. INTRODUCTION 

Magnetic field plays important roles in star formation, as well as formation and evolution of protostellar disks. 
Magnetic pressure can provide extra support to the disks and magnetic field can promote removal of angular momentum 
from disks (see Velikov 1959; Chandrasekhar 1961; Balbus & Hawley 1991). However, there are many uncertainties 
' for the structure and effects of the magnetic field in protostellar disks. 
■ ■ ' Infrared (IR) polarimetry may be an important tool to investigate magnetic field structure in protostellar disks, 
provided that the grains are aligned in the disks in respect to magnetic field. Since the grain emissivity is larger for 
the long axis of a grain, emitted radiation has a polarization vector parallel to the grain's long axis. If grains are 
aligned with their long axes perpendicular to magnetic field, the resulting grain emission has polarization directed 
perpendicular to the magnetic field. Therefore, by measuring the direction of polarization of IR emission from dust 
grains one can infer the direction of magnetic field. The key question thus is whether grain alignment is efficient in 
protostellar disks. 

The notion that the grains can be aligned in respect to magnetic field can be traced back to the discovery of star-light 
polarization by Hall (1949) and Hiltner (1949), that arises from interstellar grains. Historically the theory of the grain 
alignment was developing mostly to explain the interstellar polarization, but grain alignment is a much wider spread 
phenomenon (see Lazarian 2007 for a review) . Among the alignment mechanisms the one related to radiative torques 
(RTs) looks the most promising. We invoke it for our calculations below. 

The RTs make use of interaction of radiation with a grain to spin the grain up. The RT alignment was first discussed 
by Dolginov (1972) and Dolginov & Mytrophanov (1976). However, quantitative studies were done only in 1990's. In 
their pioneering work, Draine & Weingartner (1996, 1997, henceforth DW96 and DW97) demonstrated the efficiency 
of the RT alignment for a few arbitrary chosen irregular grains using numerical simulations. This work identified RTs 
as potentially the major agent for interstellar grain alignment. Cho & Lazarian (2005, henceforth CL05) demonstrated 
the rapid increase of radiative torque efficiency and showed that radiative alignment can naturally explain decrease of 
the degree of polarization near the centers of pre-stellar cores. Large grains are known to be present in protostellar 
disk environments and this makes the RT alignment promising. 

The effect of RTs is two- fold. They can spin- up grains and drive the alignment. While the details of the second 
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process are a subject of intensive research (see Weingartner & Draine 2003; Lazarian & Hoang 2007; Hoang & Lazarian 
2007), for our estimates we use the RTs spin-up efficiency to evaluate the efficiency of grain ahgnment. As grains of 
different temperatures are present in protosteUar disks, the differential alignment of grains at different optical depths 
is expected to show itself through variations of polarization at different wavelengths. 

ProtosteUar disks are often detected through far-infrared excess. Dust grains in the protosteUar disks are the main 
cause of the infrared excess - dust grains absorb stellar radiation and re-emit at infrared wavelengths. The spectral 
energy distribution (SED) of the emitted light gives valuable information about the disk structure. Recently proposed 
hydrostatic, radiative equilibrium passive disk model (Chiang & Goldreich 1997; Chiang et al. 2001; hereafter CG97 
and COl, respectively) fits observed SED from T Tauri stars very well and seems to be one of the most promising 
models. Here, passive disk means that active accretion effect, which might be very important in the immediate vicinity 
of the central star, is not included in the model. In this paper we use the model in COl. 

Polarization arising from disks is of great interest and importance. Recently Aitken et al. (2002) studied polarization 
that can arise from magnetized accretion disks. They considered a single grain component consisting of the 0.6/xm 
silicate and assumed that grains were partially aligned with R — 0.25, where R =< 3cos^ /3 — 1 > /2 is the Rayleigh 
reduction factor (Greenberg 1968). Here /3 is the precession angle between the grain spin axis and the magnetic field 
and the angle brackets denote average (see Aitken et al. 2002 for details). In this paper we use both a theoretically 
motivated model of grain alignment and a more sophisticated model of accretion disk. 

Note, that challenges of observing polarization in mid-IR to far-IR (FIR) wavelengths have been dealt with suc- 
cessfully recently. For instance, Tamura et al. (1999) first detected polarized emission from T Tauri stars, low mass 
protostars. As technology develops, the future of polarimetric IR and submillimeter emission studies looks very promis- 
ing. In this paper we try to theoretically predict polarized emission from T Tauri disks at different wavelengths and 
for different inclinations of the disks. 

We calculate grain alignment by radiative torque using a T Tauri disk model in COl and we predict polarized mid- 
IR/FIR/sub-millimeter emission. In §2, we discuss grain alignment in T Tauri disks. In §3, we give theoretical estimates 
for degree of polarization. In §4, we calculate the spectral energy distribution of maximally polarized emission, which 
will be useful only when we spatially resolve disks. In §5, we discuss the effect of inclination angle. In §6, we discuss 
observational implications. We give summary in §7. 

2. GRAIN ALIGNMENT IN PROTOSTELLAR DISKS 
2.1. The disk model used for this study 

We assume that magnetic field is regular and toroidal (i.e. azimuthal). We use a T Tauri disk model in COl. Figure 
[1] schematically shows the model. The disk is in hydrostatic and radiative equilibrium (see also CG97) and shows 
flaring. According to CG97, flaring of disk is essential for correct description of SED. They considered a two-layered 
disk model. Dust grains in the surface layer are heated directly by the radiation from the central star and emit their 
heat more or less isotropically. Half of the dust thermal emission immediately escapes and the other half enters into 
disk interior and heats dusts and gas there. Both CG97 and COl assumed that the disk interior is isothermal. 

In both CG97 and COl, the disk surface layer is hotter than the disk interior. Thus, roughly speaking, the surface 
layer dominates in mid-infrared wavelengths and disk interior dominates in far- infrared/sub- millimeter wavelengths. 
The disk surface layer is both optically and physically thin. 

The major difference between CG97 and COl is the treatment of dust grain size distribution. CG97 assumed that 
all grains have a fixed size of 0.1/im, while COl assumes an MRN distribution (Mathis, Rumpl, & Nordsieck 1977) 
with maximum grain size of UmaxA = lOOO/im in the disk interior and amax,s = ifJ-m in the disk surface layer. 

In our calculations, we use a grain model similar to that in COl. We use an MRN- type power- law distribution of 
grain radii a between amm (=0.01 iim for both disk interior and surface layer) and amax (=1000 fim for disk interior 
and = 1/Ltm for disk surface layer) with a power index of -3.5: dN oc a~^'^da. As in COl we assume that grain 
composition varies with distance from the central star in both disk interior and surface layer. We assume that grains 
in the surface layer are made of silicate only when the distance r is less than 6 AU, and silicate covered with water 
ice when r > 6AU. We do not use iron grains for the immediate vicinity of the star. We assume that grains in the 
disk interior are made of silicate when r < 0.8AU and ice-silicate for r > 0.8AU. The fractional thickness of the water 
ice mantle, Aa/a, is set to 0.4 for both disk surface and disk interior. Unlike COl, we use the refractive index of 
astronomical silicate (Draine & Lee 1984; Draine 1985; Loar & Draine 1993; see also Weingartner fc Dr aine 2001). We 
take optical constants of pure water ice from a NASA web site | ftp://climatel.gsfc. nasa.gov/wiscombe). 

The column density of the disk is ^arj^f/"^ with Sq — lOOOg / cm'^ . Here tau is distance measured in AU. The disk is 
geometrically flared and the height of the disk surface is set to 4 times the disk scale height h. The disk inner radius 
is 2i?* and the outer radius is lOOAU. The central star has radius of i?* = 2.5i?s„„ and temperature of T, = 4000K. 
Temperature profile, flaring of disk, and other details of the disk model are described in COl. 

2.2. Radiative torque for large grains 

For most of the ISM problems, dust grains are usually smaller than the wavelengths of interest. However, this is no 
longer true in T Tauri disks because we are dealing with grains as large as ~ lOOO/im. To understand grain alignment 
in T Tauri disks we need to know radiative torque for large grains. 

In this study, we do not directly calculate radiative torque for large grains. Instead, we use a simple scaling relation 
to model radiative torque for large grains. 
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In CL05, we used the DDSCAT software package (Draine & Flatau 1994; Draine & Flatau (astro-ph/0409262D; 
Drains & Weingartner 1996) to calculate radiative torque on grain particles and showed the relation between XQr and 
A/a for grains with radii between O.l^m and 3.2^m. Here Qr is the radiative torque efficiency. Figure^ obtained by 
reprocessing the earlier relation, shows that the radiative torque 

_r^O(l) ifA^a 
^r-|^(A/a)-3 ifA>a, 

where a is the grain size and A the wavelength of the incident radiation. Note that the radiative torque peaks near 
A ~ a and that its value is of order unity there. A more general study on this issue is provided in Lazarian & Hoang 
(2007). This allows us to assume that the relation holds true both for small and large grains. 




2.3. Rotation rate of dust grains by radiative torque 

To obtain grain rotational velocity one requires to calculate the balance between the excitation of rotation, driven 
by different processes and the damping of rotational motions (see Draine & Lazarian 1998). Below we shall calculate 
the rotational rates of a grain driven by radiative torques and subjected to gaseous damping. Note that, according to 
Lazarian & Hoang (2007), the obtained rates will not be the actual rotational rates of a grain that is free to get aligned 
under the influence of the anisotropic radiation. In fact, most of the grains will be driven to low- J attractor points and 
will rotate thermally or even sub-thermally. However, being driven by anisotropic radiation, grains get aligned and it 
was argued in Lazarian (2007) that the parameterization of the radiative torques in terms of the maximal achievable 
angular velocity is a valid one. Incidentally, this was also the parameterization that we used in CL05 and in the 
subsequent paper by Bethell et al. (2007). 




Fig. 2. — Behavior of Torque. Torque is ~ 0{l) when A ~ a, where a is the grain size. Roughly speaking, torque oc (A/a) ^. Left panel: 
The results for small grains. Data from CL05. Right panel: The results for large grains. Data from Lazarian & Hoang (2007). 
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After some modifications, equation (67) in DW96 reads 

Wr-od\^ , „„ mgO^l f '^■rad \^ f ^ \^ tr-. l2 f '''drag \^ 

= 4.72 X 10«— -p3a_5 — ; [QrV — , 2 

Wt / 0^ \nHkT J \fimj \rdrag,gas J 

where Qr = Qr • and ai is the principal axis with largest moment of inertia, nn is the hydrogen number density, 
Urad is the energy density of the radiation field, (5 « 2, ai « 1.745, ps — rho/3gcm~^, = a/10~^CTO, and lot is the 
thermal angular frequency, which is the rate at which the rotational kinetic energy of a grain is equal to kT/2. The 
timescales T^rag.gas and Tdrag,em are the damping time for gas drag and for electromagnetic emission, respectively, and 
they satisfy the relation r^^^ = '''drag em + '''drlig gas i^*^^ Draine & Weingartner (1996) for details). As we discussed 
in the previous subsection, Qr is of order of unity when A ^ a and declines as (A/a) increases. From this observation, 
we can write 

/ \ 2 / \ 2 / ^ '\ \2 / \2 

OJrad \ ^1 ^rad \ ( Vr,A~a \ ^ f l^rad \ I \ ^g-j 



,,4.72xlO''^^P3a„5 — [-^^ , (4) 

J Xr^a " \rLHkT J \fJ.mJ \Tdrag,gas J 

2.4. Minimum and maximum aligned grain size 

As in DW96, we assume that grains are aligned when (uirad/^^Jr)'^ > 10. Suppose that a monochromatic radiation 
field illuminates dust grains and that [ujrad/ > 10 for A ~ a. According to Eq. ([2]), the ratio {uJrad/^TY decrease 
as a decreases and the ratio will drop below 10. Applying Eq. ([3]) we can easily find the minimum aligned grains 
size: 

10 ^ 



UJT J V J x^a V Qr,\ J \ J 



for A > a, where 

l^rad 



9 ^1 ^ / Urad V / ( '^'^''^g 



0-lower ^ ^ 7 ; Tt I (5) 

\Kad/WT)w/ 

Then, how can we find the maximum aligned grain size? In other words, are all grains with a > A aligned? In 
order to answer this question, we need to consider the behavior of torque Qr for the limit a ^ A. Note that, when 
A ~ a, Qr ~ 0{1). Then what happens when a ^ A? In principle, we can calculate Qr using numerical simulations. 
However, this is still infeasible because enormous computational power is required. Therefore we can only conjecture 
what will happen for a ^ A. 

If one adopts the reasoning for the origin of RT in Dolginov & Mytrophanov (1976) one might expect that Qr drops 
due to incoherent contributions. Theoretical considerations in Lazarian & Hoang (2007) show that this may not be 
true. While to be conservative we adopt a rule of thumb is that Qr may begin to decline when a > lOA, we will see 
below that this assumption does not alter our results in any appreciable way^. 

We also need to consider many different time-scales: precession time-scale, time-scale for the alignment of angular 
momentum and ai, etc. However, fortunately, the exact knowledge on the maximum aligned grains size is not so 
important for our current study (see next section). 

2.5. Grain alignment in disks 

We use Eq. ([1]), instead of the DDSCAT software package, to obtain radiative torque (Qr) on grain particles in the 
T tauri disk described in ^21 We take a conservative value of Qr at A '--^ a: Qr ~ 0.1 at A ~ a. Apart from Qr, we 
also need to know Urad and uh to get the (wrad/wT)^^^ ratio (see Eq. Q). We directly calculate u^ad and nn using 
the disk model in COl. We assume that Tdrag ~ Tdrag,gas- 

In the disk interior, there are two kinds of radiation fields: one from the surface layer and the other from the disk 
interior itself. We assume that both radiation fluxes direct only along the disk vertical axis (i.e. "z" axis). As in CG97, 
we assume that half of the stellar radiation flux that reaches the disk surface enters the disk interior. The radiation 
flux from the star is ^ (a/2)(i?,/r)^fTsT^, where a is the grazing angle at which the starlight strikes the disk and i?* is 
the stellar radius, T, is the stellar temperature, and as is the Stephan-Boltzmann constant (CG97). We assume that 
the radiation flux from the surface layer has a narrow spectrum around the wavelength of ~ 3000 /T^s P-tti, where Tds 
is the grain temperature in the surface layer (see COl for a graph for Tds)- The magnitude of the flux from the surface 
at a given height is less than the half of the incident stellar flux because the flux from the surface is attenuated by 
dust absorption in disk interior. We also need to consider that there are two surface layers - one above the mid-plane 
and the other below it. The flux from disk interior is also treated as a monochromatic wave with A ^ SOOO/T^ /im. 
Note that flux from the interior has longer wavelengths because the disk interior is cooler (see COl for a graph for the 
disk interior Ti). From the two fluxes, we can calculate anisotropic component of the radiation energy density. We 

use the fact that the disk surface density is given by E = lOOOr^^^^cm"^ and that the height of the disk H{r) can be 
calculated from 



H/r ^ 4VTVT;V^7^, (6) 

^ Note that here A is the wavelength at which radiation is strongest. For most cases, A = Xmax,Wie.n 1/T when the radiation field is 
a black body radiation, where Xmax.Wien is the Xmax from Wien's law. 
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where Tc = 3 x 10^^"'GA/,//ci?* in cgs units. Here k is the Bohzmann constant. 

We assume that grains are ahgned when the ratio {tOrad/i^T)'x^a exceeds 10, which is overly conservative according 
to a recent study in Hoang & Lazarian (2007). Figure [3] shows that, at large r, large grains are aligned even deep 
inside the interior. Here we take A as the Xmax.Wien of the local blackbody radiation field at r. On the other hand, at 
small r, only grains near the disk surface are aligned. This is because gas density is low and, therefore, the gas drag is 
smaller near the surface layer. The lower panel of Figure [3] shows that almost all grains are aligned when r > lOAU. 
Therefore, we expect strong polarized emission from outer part of the disk. 

In the surface layer, grains are aligned by the star-light. Note that the radiation field scales as r~^. Since the column 
density scales as r~'^/^ and the disk height is an increasing function of r, the density in disk surface will drop faster 
than r~^/^. The gas temperature in the surface layer roughly scales as r^^"^ (see CG97 and COl). Therefore we expect 
that the ratio oJrad/i^T is an increasing function of r (see Eq. ([2|)). This implies that grains at large r are aligned. 
Grains near the central star cannot be aligned due to high gas density near the star. Indeed Figure H] shows the ratio 
{'^rad/i^T)'\^a exceeds 10 when r > lAU, which means that grains in the surface layer are aligned when r > lAU. 
Note that the radiation from the central star has a Xmax,Wien at ~ 750nm. We expect that polarized emission from 
the surface layer is also originated from outer part of the disk. 




0.01 0.10 1.00 10.00 100.00 0.01 0.10 1.00 10.00 100.00 

log r(AU) log r(AU) 

Fig. 3. — Grain alignment in disk interior. Upper panel: Contours show the ratio {u)^ad/^T)\^a^ '^o*^ that the disk vertical height is 
shown in units of the disk scale height h. We assume that grains are aligned when {Lj^^d/'^T)'^ > 10- Lower panel: Fraction of grains that 
satisfy {wj.ad/'-^T)\r^a > 10 as a function of disk radius. After r > lOAU, almost 100% of grains satisfy the inequality. 

Fig. 4. — Grain alignment in surface layer. The ratio (uj^^j^/ ujx)\^a exceeds 10 when r > lAU, which means that some grains in the 
surface layer are aligned when r > lAU. 



3. THEORETICAL ESTIMATES OF DEGREE OF POLARIZATION 

In this section, we estimate the degree of polarization of emitted radiation in IR wavelengths. We need only the 
grain size distribution for the calculation in this section. We do not use a detailed disk model. The discussion in this 
section is applicable for any system with large grains. 

Suppose that all grains are perfectly aligned. Then, what will be the observed degree of polarization? Of course, we 
do not observe 100% polarization. The size parameter 27ra/A plays an important role here. There are three important 
points that determine the degree of polarization of emitted radiation. 

First, when the grain is small compared with the wavelength (i.e. when 27ra/A < 1), grain's intrinsic shape gives 
a limit. Suppose that grains are oblate spheroid, that the symmetry axis of the grain is parallel to y-axis, and that 
radiation is propagating along x-axis. When we send two linearly polarized radiation fields with electric field parallel 
to and perpendicular to the grain's symmetry axis, respectively, the radiation fields experience different cross-sections: 
the radiation with E || ai sees smaller cross-section. Here E is the electric field of the radiation field and ai is the 
grain's symmetry axis. As a result, we observe polarization because two cross-sections are different. So far we have 
dealt with polarization by absorption. Polarization by emission is caused by the exactly same fact that two cross- 
sections are different. Since Qabs = Qerm where Qem is the grain emissivity, grain emits more radiation with E _L ai 
than one with E || ai The degree of polarization of emitted light is {Qabs,i. — Qabs,\\) I {Qabs.i^ + Qabs,\\)j where || and 
_L refer to directions parallel and perpendicular to the grain's symmetry axis ai. 

Figure El obtained from the DDSCAT package, shows this effect clearly. The the radiation fields have A = 850/im 
and the grains' long to short axis ratio is 1.3:1. The ratio of two cross-sections is around 1.6 for 27ra/A < 1. When 
we observe emission from those grains, the degree of polarization can be as large as (1.6-1)/(1.6-|-1)=22%. If we 
assume that grains' long to short axis ratio is 1.5:1, the ratio of the cross-sections is 2.1:1 and the resulting degree of 
polarization for emission is as large as (2.1-1)/(2.1+1)=35%. The ratio of the cross-sections varies as wavelength of 
radiation varies. The ratio seems to be fairly constant for A > lOO/im. However, for shorter wavelengths, the ratio 
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0.1 1.0 10.0 100.0 1000.0 0.1 1.0 10.0 100.0 1000.0 

Grain Size (/.im) Grain Size (/.im) 

Fig. 5. — Relation between the absorption cross-section {=Qabs times the geometric cross-section) and grain size. The wavelength at 
which we observe is 850 fim. The grains are oblate spheroid made of silicate with the symmetry axis ai parallel to y-axis. Axis ratio is 
1:1.3. Radiation is propagating along x-axis. Left panel: When grains are smaller than ~ 850/27r, then we can observe polarization because 
two cross-sections are different. When grains are larger than ~ 850/27r, then we cannot observe polarization because two cross-sections are 
similar. Right panel: When the axis ratio is 1:1.3, the ratio of Qabs is around ~ 1.6 when the grain is small compared with the wavelength 
(i.e. 27ra/A < 1) and ~ 1 for geometrical optics regime (i.e. 27ra/A > 1). 
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Fig. 6. — The ratio Qabs,±/Qab3,\\ versus wavelength. Left panel: Silicate grains. For A > lOO/im, the ratio is ~ 1.6 when the size 
parameter is less than 1 (i.e. 27ra/A < 1). However, for shorter wavelengths, the ratio for 2iTa/X < 1 drops: it is ~ 1.3 for A = lOfim. Right 
panel: Ice-Silicate grains (i.e. grains with silicate core and water ice mantle). 

decreases (Figure [6]). The shape of grains in protosteUar disk is uncertain (see, for example, Hildebrand & Dragovan 
1995 for the general ISM cloud cases). 

Second, when 27ra/A > 1 (i.e. in the geometrical optics regime), we do not observe polarization. Figure O clearly 
shows that the ratio of cross-sections becomes very close to 1 when 27ra/A > 1. This means that the usual argument 
about polarization by absorption or emission works only when the grain size is small compared with the wavelength: 
27ra/A < 1. That is, in the small size parameter case, when radiation meets an elongated grain, it recognizes the 
elongated shape and interacts differently depending on the direction of the electric field of the radiation. However, in 
the geometrical optics regime (i.e. when 27ra/A > 1), radiation does not recognize that grains are elongated. We can 
easily understand this fact when we consider an elongated macroscopic object: Cross-sections are same regardless of 
the electric field directions for visible light. This second point is somewhat tricky: Even in the case grains are actually 
"aligned" , we do not "observe" polarization when the grains are large compared with the wavelength. 

This observation has an important consequence. When we calculate polarization, large grains (i.e. grains with 
a > A/27r) do not contribute to polarization even when they are aligned. This fact reduces the degree of polarization 
significantly in protosteUar disks. For example, suppose that we has grains as large as lOOO/im in a disk. If we observe 
emission from the disk at A = 850/iTO, grains with a > 850/27r ~ 100/im do not contribute polarization although they 
dominate extinction when the grain size distribution is a power-law (dN oc a~'^da) with g < 3. 

Thirds, if the medium is optically thick, the degree of polarization reduces. The intensity of radiation from a uniform 
slab is Sv{l — ea;p(— r)), where Si, is the source function and r is the optical depth. We can observe polarization when 
r, is different for parallel and perpendicular (to grain symmetry axis or generally any spatial direction) directions. 
However, in the optically thick limit, the intensity becomes equal to 5*,^. Therefore, difference in optical depth does 
not produce observable level of polarization when the slab is opaque. 
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For the disk interior^ let us consider only the first and the second points mentioned above. That is, let us assume 
that the grains are perfectly aligned and the disk is optically thin. We can estimate the degree of polarization p for 
the disk interior by integrating the following: 

^^^^ _ lai^^^ {QabsM - Qab.,\\{a)) a^N{a)da 
XTr.r (<3abs,_L(a) + Qabs,\\{a)) a'^N{a)da 

where Q... is grain absorption efficiency, or grain emissivity, and amin — 0.01/im. In left panel of Figure [71 we use 
O'max = 1000/iTO. Wc usc dN oc a~'^da with q — 3.5 and assume that grains are oblate spheroid with long to short 
axis ratio of 1.5:1. It is not surprising that the degree of polarization rises when we use smaller amax- when amax is 
smaller, more grains are in the geometrical optics regime. Right panel of Figure [7| shows this effect. 

Note that we do not use actual disk models here. Actual numerical simulations using actual disk models will give 
smaller values because not all grains are aligned and some part of the disk is optically thick. 

Left panel of Figure [7] shows that the maximum degree of polarization for A < lOO/im is slightly larger than that 
for A > 100/im. However, in reality, we do not expect a significant degree of polarization for A < 100/im. The reason 
is that the entire disk becomes optically thick for A < 100/im. The opacity per unit mass for A = 100/im is around 
~ 0(0. 1) (see Figure 3 in COl). The outer most disk has column density of Ig/cm^. Therefore, even the outer-most 
part of the disk becomes optically thick when the wavelength drops below ~ 100/im. 

For the disk surface layer, the second point mentioned above is irrelevant because grains are smaller in the surface 
layer: amax — IfJ-m. The third point is also irrelevant because the surface layer is optically thin at far-infrared and 
submillimeter wavelengths. Therefore, if all grains are perfectly aligned, the degree of polarization of the emitted 
radiation is determined only by the grain shapes (see the first point above). 




10 100 1000 200 400 600 800 1000 

Wavelength (/xm) Maximum Grain Size a^^^ if^^) 

Fig. 7. — Expected maximum degree of polarization for disk interior. Left panel: amax = lOOOfim. We assume grains are oblate 
spheroid with axis ratio of 1.5:1. Right panel: The degree of polarization is very sensitive to the maximum grain size, amax- Results are 
for ice-silicate. When amax is smaller, more grains are in the geometrical optics regime and, hence, the degree of polarization rises for a 
given observing wavelength. 



4. ESTIMATES FOR SPECTRAL ENERGY DISTRIBUTION 

In this section, we calculate the degree of polarization of emitted infrared radiation from a disk with structure and 
parameters described in COl. In this section, we assume that the disk in face-on. The degree of polarization will 
be zero for a face-on disk when magnetic field is perfectly azimuthal and the disk is cylindrically symmetric. In this 
section, we are concerned only with the absolute magnitude of the polarization. 

4.1. Spectral energy distribution 
When (uJrad/^T)'\r^a larger than 10, we find the minimum aligned size from Eq. ([5]): 

where Xmax,Wien IS the peak wavelength of the aligning radiation. Grains smaller than a/otoer are not aligned. 
When (wr.ad/i^T)A^a IS larger than 10, we find the maximum size of grains that give rise to polarization from 

^upper — ^^^[10A77iaa;,Wie?i; Ao5s/27r], (9) 

where Xobs the observing wavelength. In most cases, aupper — Aofcs/27r because Xmax,Wien falls in far-infrared wave- 
lengths. Note again that the actual maximum aligned size can be larger than aupper- 



(Slower ^ \ ( , I . . \2 ) '^maK.Wiem (8) 
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Now we know a^. 



and aio 



Note that a^. 



and aio 



are functions of the distance to the central star, r, 



and the distance to the disk mid-plane, z. We calculate the parallel (with respect to the local magnetic field) and 
perpendicular opacity respectively: 

(10) 

(11) 



T|[(r,z)cx / Qabs{T^a^)N{a)f\\da, 

Tj^(r, z)cx / Qabs{'^a'^)N{a)f±da, 

where N{a)da cx a~^-^ (i.e. an MRN-type distribution) and 

J. — ^''^'^ ^ 0.65) if aiower < a < ttu 

and 



otherwise, 



„ _ J - 1.23 (or - 1.35) if 
■'^ 1 1 otherwise, 

where we assume that the long to short axis ratio of the oblate spheroid is 1.3:1 (or 1.5:1). 
From this, we can calculate emission in parallel and perpendicular directions: 



oc A 



L\,j_ oc A 



ih 



dT 



A, II 



dr r 1 

J-4h 

-4h 



dz^e--^Bx{T), 



(12) 



(13) 



(14) 



(15) 



where r|| and T|[ measure optical depths from z to Ah along the axis perpendicular to the disk mid-plane (see CG97). 
We use BHCOAT.f and BHMIE.f codes in Bohren & Huffman (1983) to calculate grain emissivity Qabs- 
The degree of polarization is 

p{\) = (iAa - iA,||)/(iA,± + iA,||). (16) 

Figure [5] ( Ze/i? panel for 1.3:1 oblate spheroid and right panel for 1.5:1 oblate spheroid) shows the results. The degree of 
polarization can be as large as ~5% in FIR/sub-millimeter wavelengths and 10% in mid-IR regimes. The polarized 
emission at FIR is dominated by the disk interior and that at mid-IR is dominated by the disk surface layer. Note 
again that, in these calculations, we ignored the direction of polarization and we only take the absolute value of it. 

4.2. Radial energy distribution 

Figure O shows radial distribution of emitted radiation. For A = 850/^m, both radiations from the disk interior and 
the surface layer are dominated by the outer part of the disk. But, for A = 10/im, the inner part of the disk contributes 
significant portion of total emission and polarized emission. 

5. EFFECTS OF DISK INCLINATION 

In this section we calculate actual degree of polarization that we can observe. Chiang & Goldreich (1999) calculated 
spectral energy distribution (SED) from inclined disks. We follow a similar method to calculate the the SED of 
polarized emission. The SED of disk interior is the integral of the source function: 



A 



dx 



dy / dTxBx{T,)e' 



where 



y 



i0 - H{rrnax) sin 0, 



(17) 



(18) 



and rmax is the outer disk radius, H{r„iax) is the height of disk at the disk outer radius, is the angle between disk 
symmetry axis and the line of sight, Ti is the temperature of disk interior, Bx is the Planck function, and t is the 
optical depth. The integral / drx--- is taken over the line of sight (see Chiang & Goldreich 1999 for details). The SED 
of the disk surface is obtained from the integral 



L 



surf 



OC A 



dx / dyy^BxjTds 



y(x) 



1 — exp 



In- 11 



exp{-Tx), 



(19) 



where fi and 1 are unit vectors normal to the surface and parallel to the line of sight, respectively, Cg is the Planck 
averaged dust emissivity at the surface. The summation is performed whenever the line of sight intersects the surface 
(see Chiang & Goldreich 1999 for details). 

In our calculations, we explicitly take care of the fact that grain symmetric axis is changing along a given line of 
sight. We follow the description in Roberge & Lazarian (1999) (see also Lee & Draine, 1985) to calculate this effect 
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Fig. 8. — Spectral energy distribution. The vertical axis (i.e. A_Fx) is in arbitrary unit. Thick solid line: total (i.e. interior + surface) 
emission from disk. Thin dotted line: total emission from disk surface. Thick dotted line: polarized emission from disk surface. Thin 
dashed line: total emission from disk interior. Thick dashed line: polarized emission from disk interior. Note that, in these calculations 
of polarized emission, we ignored the direction of polarization vectors and we only take the absolute value of them. Left panel: Results 
for oblate spheroid grains with axis ratio of 1.3:1. Right panel: Results for oblate spheroid grains with axis ratio of 1.5:1. The degree of 
polarization is larger than that in the left panel. 

and we obtain optical depths with respect to the x and y directions in the above integrals (Eqs. pT|) and (|19p )^. After 
calculating the optical depths, we calculate Lx ^ and L\^y. We obtain the luminosity for two more directions, which 
have 45 degrees with respect to x- and y-directions. We calculate the the direction and degree of polarization based 
on the luminosity for these 4 directions. 

Figure [10] shows the effects of the disk inclination. We calculate the polarized emission from the disk interior. 
The viewing angle 6 (=the angle of disk inclination) is the angle between the disk symmetry axis and the line of 
sight. We plot the direction of polarization for 3 different wavelengths and 2 different viewing angles. The lines 
represent the direction of polarization. Since we assume that magnetic field is azimuthal, the direction of polarization 
is predominantly radial (see lower panels). In Figure [TT] we show similar plots for radiation from the disk interior only. 
For A > 100/im, the polarization patterns in Figure [TT] is very similar to those in Figure [TOl But near the disk edges. 
Figure [To] shows larger degree of polarization than Figure [TT] This is because the emission from the disk interior is very 
weak there compared with that from the disk surface layer. For A < 100/j,to, the polarization patterns in Figure [TT] 
are very different from those in Figure [TOl because polarized emission from the disk surface layer dominates that from 
the disk interior. Note that, since the degree of polarization of emission from the disk surface layer is very sensitive 
to the maximum grain size in the surface layer, the results for A < lOO/im should be very sensitive to the maximum 
grain size in the surface layer. 

While the polarimetry of the spatially resolved accretion disks is promising with new generation of instruments (see 
§6.2), at present one can study disk magnetic fields with unresolved accretion disks. Below we provide predictions for 
this case. Figure [T2l shows spectral energy distribution for such a disk for four different viewing angles. When 9 = 90 

^ Let the z-axis be the direction of the line of sight and the yo axis the direction along the projection of the magnetic field on to the 
plane of the sky. The xq axis is perpendicular to both axes. Then, the cross sections CxO a-nd Cyo are 

c,o = Cavg + lRic± - c-y)!! - SCOS^C), 

CyO = Cavg + —R{C± — C||), 

where Cavg = (2Cx + C||)/3, C± and C|| are cross sections with respect to the magnetic field, is the angle between the magnetic field 
and the plane of the sky. We assume the Rayleigh reduction factor, R, is 1. For the axes x and y that also lie in the plane of the sky, 

Cx ~ CxO cos^ 6 + CyO sin^ 9, 
Cy Rj CxO sin^ 9 + CyO cos^ 9, 

where 9 is the angle between the x and xo axis. 
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Fig. 9. — Radial energy distribution, (a) A = lO/im. Inner part of the disk emits substantial amount of radiation. But it emits negligible 
amount of polarized radiation. Note that, when r < lAU, grains in the surface layer are not aligned and only negligible fraction of grains 
are aligned in the interior (see Figures |3] and |4]l . (b) A = bOfim. (c) A = lOO/^m. (d) A = 850/jm. The result for A = 450/xm (not shown) is 
very similar to that for A = 850/nn. 

(i.e. for edge-on disk), inner part of the disk (i.e. region close to the star) is invisible due to high opacity. Therefore 
the spectral energy distribution truncates for A < lO/im. When 9 = (i.e. for face-on disk), the polarized emission is 
zero as expected'^. 

Finally, Figure [T3l shows the change of the degree of polarization for selected wavelengths. Left panel shows the 
degree of polarization for total emission, while the right panel shows that for radiation from the interior only. The 
degree of polarization is large when the angle is small (see left panel). The sudden drop for A = 10 fim is due to 
the following reason. As the viewing angle drops, the inner part of the disk suddenly becomes visible, which causes a 
sudden increase of the the total intensity (or flux). But, the polarized intensity (or flux) does not change much because 
the inner part of disk does not emit polarized emission. Note that grains are not aligned in the inner part of the disk. 
For A = 5G/im, polarization is dominated by the disk surface layers. 

6. DISCUSSION 
6.1. Other Alignment Processes 

In the paper above we concentrated on only one alignment process, namely, on the RT mechanism. This mechanism 
has became so promising in explaining of interstellar polarization partially because the main competitor, namely, 
paramagnetic alignment mechanism (Davis-Greenstein 1951) and its later modifications were shown to have problems 
with aligning interstellar grains. In fact, the fast spin-up mostly due to H2 formation on the catalytic sites over 
grain surface as suggested in Purcell (1979) is a textbook process that is invoked to explain the efficient paramagnetic 
alignment. Indeed, fast rotating grains should be immune to randomization by atomic bombardment and thus get 
aligned well. However, the Purcell's spin-up was shown to be inefficient for most of grains in diffuse interstellar gas due 
to thermal flipping of grains that was reported in Lazarian & Draine (1999ab). The thermal flips arise from the coupling 
of vibrational and rotational degrees of freedom stemming from the processes of internal relaxation within interstellar 
grains, in particular, Barnett (Purcell 1979) and nuclear relaxation (Lazarian & Draine 1999b). As a result of flipping 
the direction of the Purcell's torques acting on a grain alters and the grain gets "thermally trapped" . It rotates at a 
thermal velocity and therefore is subjected the randomization due to random gaseous and ionic bombardment. 

The difference of the interstellar and disk grains is their size. The larger grains in the disk are not thermally trapped. 
Therefore, potentially, the processes of Purcell's spin-up are applicable. If temperatures of grains and ambient gas are 
different, this may result in a spin- up that arises from variations of the accommodation coefficient over grain surface 

^ That is, we do not see thick dotted or thick dashed lines in Figure [T2l d) . 
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X = 850 iJ.m X=100 fim X=10 fim 




0=10° 10% 0=10" 10% 0=10° 10% 



Fig. 10. — Simulated observations. Degree of polarization is calculated for the total radiation (i.e. interior + surface) from the disk. 
The disk inclination angle 9 is the angle between disk symmetry axis and the line of sight. 

(Purcell 1979), provided that the temperatures of the gas and the grains in the disk differ. The problem of such a 
scenario is that the paramagnetic alignment is slow, unless grains demonstrate enhanced magnetic susceptibility (e.g. 
arc super-paramagnetic) (sec Jones & Spitzcr 1967). Potentially, this process that also aligns grains with long axes 
perpendicular to magnetic field can enhance the alignment and therefore polarization. However, we do not know about 
the abundance of the required grains within the disks. 

For the largest grains, a particular mechanical alignment, which was termed in Lazarian (1994) "weathercock mech- 
anism" is applicable. In the presence of gas-grain motions large irregular grains would tend to get aligned with long 
dimension along the flow as their center of pressure and center of mass do not coincide. However, the mechanism 
requires substantial relative velocities of gas and grain, which is not certain in the protostellar disks. Moreover, we 
have showed above that very large grains do not produce polarized radiation at the wavelength we deal with in this 
paper. 

More promising may be the alignment of helical grains first mentioned in Lazarian (1995) and discussed at some depth 
in Lazarian & Hoang (2007). The mechanism is based on the interaction of an irregular grain with a flow of atoms. 
If some fraction of colliding atoms is not being absorbed by the grain surface due to grain accommodation coefhcient 
not being equal to unity, the collisions with a flow of gaseous atoms should cause the alignment similar to that by 
anisotropic radiative torques. However, as it is discussed in Lazarian (2007) we do not any compelling evidence for the 
operation of the mechanism in the studied astrophysical environments. On the contrary, the alignment by radiative 
has been proven to provide the observed alignment in a number of circumstances (see Lazarian 2007). Thus we defer a 
quantitative discussion of the mechanical alignment of helical grains. If the mechanism operates efficiently, it can only 
increase the degree of alignment at the parts of the disk where the radiative torques start to fail, making polarization 
from aligned grains only more important. 

All in all, while further studies of alternative alignment mechanisms seems necessary, at present the discussed RT 
mechanism provides the safest bet. 

6.2. Observational Prospects 

Multifrcquency observations of protostellar disks have become a booming field recently. They have advanced sub- 
stantially our knowledge of the disks and allowed theoretical expectations to be tested. 

Magnetic fields are an essential components of the protostellar disks. They are likely to be responsible to accretion 
(see Nomura 2002). Therefore observational studies of them are essential. In this respect our paper is the first, as far 
as we know, attempt to provide the expectations of the polarization arising from accretion disks that is based on the 
predictions of the grain alignment theory. 
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X = 850 fj.m X=100 /j.m X=10 /^m 




0=10° 10% 0=10° 10% 0=10° 10% 



Fig. 11. — Simulated observations. Degree of polarization is calculated for the radiation from the disk interior only. 



Our study reveals that multifrequency polarimetry is very important for the protostellar disks. The synthetic 
observations that we provide expUcitly show that observations at wavelength less than 100 /im mostly test magnetic 
fields of the skin layers, while at longer wavelengths test magnetic fields of the bulk of the disk. Therefore polarimetry 
can, for instance, test theories of accretion, e.g. layered accretion (Gammie 1996). 

Combining the far-infrared polarimetry with polarimetric measurements at different frequencies may provide addi- 
tional insight into the magnetic properties of protostellar accretion disks. For instance, circular polarization arising 
from scattering of starlight from aligned grains (see Lazarian 2003) and polarization in emission lines arising from the 
aligned atoms (see Yan & Lazarian 2006, 2007) can provide additional information about the magnetic in the outer 
parts and above the accretions disks. 

Most of the present day polarimetry will be done for not resolved protostellar disks. The size of the T Tauri disks is 
usually less than ~ 300 AU (see, for example, COl). If we take the distance to proto-stars to be around > lOOpc, then 
the angular sizes of the disks are usually smaller than 6". The angular resolution of SCUBA polarimeter (SCUPOL) is 
around 14" (Greaves et al. 2000) and that of SHARC II polarimeter (SHARP; Novak et al. 2004) at 350^m is around 
9". Therefore it is not easy to obtain plots like Figures fTOl The angular resolution of the intended SOFIA polarimeter 
is around 5" at 53^m, 9" at 88^m, and 22" at 215^m. We see that the intended SOFIA polarimeter will be at the 
edge of resolving structure of close-by disks, while other instruments will not resolve typical T Tauri disk. Therefore 
for most of the near future observations our predictions in Fig. 13 and 14 are most relevant. 

Higher resolution polarimetry is expected in future, however. This will make our predictions of polarization the 
resolved accretion disks in Fig. 11 and 12 testable. Note, that the actual structure of magnetic fields may be much 
more complex than the one in our simple model. 

While in the paper we dealt with the protostellar accretion disks, our results are suggestive of importance of 
polarimetric studies of magnetic fields in the disks of evolved stars. In a broader context, the present paper is one of 
the first studies that make use of the advances in grain alignment theory to extend the utility of polarimetric studies 
of magnetic field beyond its traditional interstellar domain. 

6.3. Effects of scattering 

In this paper, we do not consider the effects of scattering. Below we show that effects of scattering is indeed less 
important for emissions from disk interior in the FIR regime. We will provide more detailed calculations on the effects 
of scattering for shorter wavelengths elsewhere. 

In the case l-najX ^ 1 we can assume that Qscatt is negligibly smaller than Qabs- However, since we are dealing 
with large grains here, we cannot simply assume the inequality. We first compare the relative magnitudes of mass 
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Fig. 12. — Spectral energy distribution for four different viewing angles. When 8 = 90 (i.e. for edge-on disk), inner part of the disk (i.e. 
region close to the star) is invisible because it is occulted by the outer part of the disk. Therefore the spectral energy distribution truncates 
for A < lOfim. When 9 = (i.e. for face-on disk), the polarized emission is zero as expected. This is because the assumed magnetic field 
configuration is perfectly azimuthal. 
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Fig. 13. — Degree of polarization vs. viewing angle. Left panel shows the degree of polarization for total (i.e. interior + surface) 
emission, while the right panel shows that for radiation from disk interior only. 

scattering and absorption coefficients, Kscatt and Kabs, when grain size distribution follows an MRN-type power law, 
dN oc a^^^^da, between amin = 0.01/im and Qmax (see also section 2.1). Figure [T4l shows Kscatt is larger than Kabs for 
A > 100/im in case amax — lOOO/im. This is because the scattering efficiency is larger than the absorption one when 
A ~ 27ra (Fig. I15[). This trend is generally observed when A > 100//m. When the cutoff size, Umax, gets smaller, Kscatt 
becomes subdominant (see the right panel of Fig. [T4|) . 

Polarization by scattering is proportional to ~ J\Kscatt, where J\ is the mean radiation field, while that by emission 
is proportional BxKabs, where Bx is the intensity of the blackbody radiation at the point of interest. Figure [TBI 
shows this J Kscatt to BKabs ratio at selected points on the disk mid-plane. We use the disk model in COl to calculate 
J\ and B\. We only include emission from disk interior. Recall that most FIR emission from disk interior are from 
r > lOAU (see Fig. IH]). Therefore, the J Kscatt to BKabs ratio for r > lOAU concerns us most. This crude estimate tells 
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Fig. 14. — Mass absorption and scattering coefficients. Y-axes are in arbitrary units. Left panel: The maximum grain size, amax, is 
1000/^m. The scattering coefficient is larger when A > 100/jm. Right panel: The maximum grain size, amax, is 50/im. The scattering 
coefficient is similar to or less than the absorption coefficient. In both cases, the minimum grain size is O.Olfim. 

US that polarization by scattering is less important than that by emission in the FIR regime. Nevertheless, readers 
should keep in mind that our polarization calculations in this paper reflect only the emitted component and, therefore, 
the radiative transfer adopted is approximate. 
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Fig. 15. — Scattering and absorption coefficients. The wavelength is fixed for A = 850/xm. The grains are spherical silicate coated with 
water ice. The scattering coefficient is larger than the absorption one for large grains. 

Fig. 16. — Comparison between polarization by scattering and emission. When the value in y-axis is less than 1, we expect that absorption 
is more important for polarization. The values are calculated in the mid-plane of the disk. When we calculate the values at locations above 
the mid-plane the values may be smaller. 



How important is the self scattering effect? Let us first consider the case that the entire disk volume has aligned 
grains (Case I in Fig. [T7|) . As we will see right below, self scattering/absorption does not change the direction of 
polarization in Case I in Fig. [T71 The observed radiation parallel to the magnetic field becomes 



-^A,!! = J Bxexp{-T^^t,\\)dTabs,\\ 



. Bxcxp{-Text,\\)d{Text.\\) ^ — [1 - exp(-re:rt.||)], (20) 

1 + a|| J 1 + Q!|| 

where Tf.^t,\\ = ''afcs, || + '''scatt.w and we assume that the gas is uniform and that a|| = Tscatt, || /Tabs. || does not change 
along the line of sight. For simplicity, we ignore the radiation that scatters into the line of sight. We will have a similar 
expression for the perpendicular direction. Therefore, we have 

1 - exp(-Te^ta) 1 -cxp(-re^tji) 

^ ^11 °^ TT ' TT ' ^^'^f 

1 -I- a± 1 + ail 

which is positive when Text.± > Te^:*,!! and a± ~ a||. This means that the direction of polarization is perpendicular 
to the magnetic field. Recall that most of the outer part of disk (i.e. r > lOAU) has aligned grains (see Figure [3]). 
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Therefore, the situation is similar to Case I and we expect that including self scattering does not change our qualitative 
results in earlier sections. 

In Case II, however, self scattering can be potentially important. Case II will happen when we have a slab of aligned 
grains in front of unaligned grains (see Fig. I17p . The observed intensity becomes 



Ix,\\ = I\,bg exp(-re^tji) + YTpT^t^ ~ exp(-re^tjl)] 
for the direction parallel to the magnetic field, and 



h,± = h.bg cxp(-re, 



+ 



1 + a_L 



[1 - CX.p{-Text,±)] 



(22) 



(23) 



for the direction perpendicular to the magnetic field. Here, for simplicity, we assume the medium is uniform. In 
optically thin case, we have 



B 



A 



1 + a 



-Text,\\ — I\,bg ~ h,bgText,\\ + B\Tabs,\\ 



— h,bg - IxMgText, 

for the direction parallel to the magnetic field, and 

Bx 

Ix,± = Ix.bg — Ix,bg^ext,± + ^ _^ ^ T^xt,!^ = Ix,bg — Ix,bgText,± + B\Tabs,± 

for the direction perpendicular to the magnetic field. Then we have 

- /|| CX Bx{Tabs,± - Tabs,\\) - Ix,bg{Text.± - Text,\\) 
'^scatt^l^ '^scatt, 



(24) 



(25) 



Bx - hj 



{Tabs,± - TabsAl)- 



(26) 



Tabs,± - Tabs,\\ 

Therefore, the direction of polarization can change as a result of scattering. 

The ratio {Tscatt,± — Tscatt.\\)/i'''abs,± — '''abs,\\) caunot be very large. Fig. \TE\ implies that scattering does not cause 
polarization in the geometrical optics regime (i.e. X/2tt < a). On the other hand. Fig. [T51 shows that scattering is 
negligible when a is a few time smaller than (A/27r). This means that scattering dominates polarization only for a 
hmited range: /3(A/27r) < a < (A/27r), where /3 w 0.2 — 0.3. However, scattering dominates absorption for /3(A/27r) < a 
(see Fig. [T5|) . Therefore, we expect that 

(27) 



'^scatt,± '^scattAl '^scatt 
■ < 



'^abs,!. '^abs,\\ 



'^abs 



(see Fig. [Hfor the ratio Kscatt/i^ab 



B 



iv,bg- 



B 



Case I 



Case II 



Fig. 17. — In Case I, all parts of the disk along the line of sight have aligned grains. In this case, the direction of polarization is 
perpendicular to the magnetic field. In Case II, the aligned part lies in front of unaligned part. In this case, the direction of polarization 
depends on the relative strength of the local source function and the background radiation. Note that I^fig is unpolarized radiation from 
background regions of the disk. 



7. SUMMARY 

Making use of the recent advances in grain alignment theory we calculated grain alignment by RTs in a magnetized 
T Tauri disk. Based on this, we calculated polarized emission from the disk. Our results show that 

• Polarization arising from aligned grains reveals magnetic fields of the T Tauri disk. 

• Grain size distribution is the most important factor that determine the degree of polarization. 
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Fig. 18. — Scattering efficiency vs. grain size. When grains are oblate spheroid, the scattering cross-section depends on the direction of 
polarization, similar to that of absorption (see Fig. 5). See the caption of Figure 5 for details. 

• Disk interior dominates polarized emission in FIR/sub- millimeter wavelengths. When there are many grains with 
maximum grain size of lOOO/xm, the degree of polarization is around or less than ^ 2% in these wavelengths. 
However, when the maximum grains are smaller, we expect higher degree of polarization. 

• Disk surface layer dominates polarized emission in mid-IR wavelengths. The degree of polarization is very 
sensitive to the maximum size of grain in the disk surface layer. When the maximum grain size is as large as 
~ Ifim, we expect ~ 10% of polarization at A ~ 50/xm. However, when the maximum size is smaller the the 
value will drop. 

• Our study of the effect of the disk inclination predicts substantial changes of the degree of polarization with 
the viewing angle. The coming mid-lR/FIR polarimeters are very promising for studies of magnetic fields in 
protostellar disks. 
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